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Abstract 

The Langevin system subjected to non-Gaussian noise has been discussed, by 
using the second-order moment approach with two kinds of models for generating 
the noise. We have derived the effective differential equation (DE) for a variable x, 
from which the stationary probability distribution P{x) has been calculated with 
the use of the Fokker-Planck equation. The result of P{x) calculated by the moment 
method is compared to several expressions obtained by different methods such as 
the universal colored noise approximation (UCNA) [Jung and Hanggi, Phys. Rev. 
A 35, 4464 (1987)] and the functional-integral method. It has been shown that 
our P(x) is in good agreement with that of direct simulations (DSs). We have also 
discussed dynamical properties of the model with an external input, solving DEs in 
the moment method. 
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1 INTRODUCTION 



Interesting, unconventional phenomena such as the stochastic resonance (SR) and the 
noise-induced phase transition are created by noise. Theoretical studies on noise in non- 
linear dynamical systems have usually adopted Gaussian white (or colored) noise. In 
recent years, there is a growing interest in studying dynamical systems driven by non- 
Gaussian noise. This is motivated by the fact that non- Gaussian noise with random 
amplitudes following the power-law distribution is quite ubiquitous in natural phenom- 
ena. For example, experimental results for crayfish and rat skin offer strong indication 
that there could be non-Gaussian noise in these sensory systems [JJ[2]. A simple mecha- 
nism has been proposed to generate the non-Gaussian noise [3]. With the use of such a 
theoretical model, the SR induced by non-Gaussian colored noise has been investigated 
[I]. It has been shown that the peak in the signal-to- noise ratio (SNR) for non-Gaussian 
noise becomes broader than that for Gaussian noise. This result has been confirmed by 
an analog experiment [5]. 

Stochastic systems with non-Gaussian colored noise are originally expressed by the 
non-Markovian process. This problem is transformed into a Markovian one by extend- 
ing the number of variables and equations. The relevant Fokker-Planck equation (FPE) 
includes the probability distribution expressed in terms of multi-variables. We may trans- 
form this FPE for multivariate probability to the effective single-variable FPE, or ob- 
tain one-variable differential equation (DE) with the use of some approximation methods 
like the universal colored noise approximation (UCNA) [6l [7] and the functional-integral 
methods [8] [9] . The obtained results, however, do not agree each other, depending on the 
adopted approximations, as will be explained in Sec. 2.2 (Table 1). It is not easy to trace 
the origin of this discrepancy because of the complexity in adopted procedures. The pur- 
pose of the present paper is to discuss the non-Gaussian noise and to make a comparison 
among various methods, by employing the second-order moment method which is simple 
and transparent, and which is exact in the weak-noise limit. 

The paper is organized as follows. We have applied the second-moment method to 
the Langevin model subjected to non-Gaussian noise which is generated by two kinds 
of models. In Sec. 2, non-Gaussian noise is generated by the specific function which 
was proposed by Borland [3] and which has been adopted in several studies [4][8J|J]. In 
contrast, in Sec. 3, non-Gaussian noise is generated by multiplicative noise [rO]-[H]. We 
derive the effective one-variable DE, from which the stationary distribution is calculated 
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with the use of the FPE. A comparison among various methods generating the non- 
Gaussian noise is made is Sec. 4, where contributions from higher moments than the 
second moment are also discussed. The final Sec. 5 is devoted to our conclusion. 



2 Models A and A 
2.1 Moment method 

We have adopted the Langevin model subjected to non-Gaussian colored noise (e) and 
Gaussian white noise (V>£)> as gi ven by [3] 

x = F (x) + e(t)+W(t) + I(t), (1) 
re = K(e)+(jyr](t), (model A ) (2) 

with 

m = "[l + (g-l)(r/W (3) 
which is referred to as the model Ao- In Eqs. (l)-(3), F(x) is an arbitrary function of 

x, I(t) stands for an external input, q is a parameter expressing a departure from the 

Gaussian distribution which is realized for q = 1, r denotes the characteristic time of 

colored noise, and rj and £ the zero-mean white noises with correlations: (f](t)i](t')) = 

6(t - 0, mm) = S(t - and {n(t)ffl) = 0. 

First, we briefly discuss the non-Gaussian colored noise generated by Eqs. (2) and (3), 

which yield the stationary distribution given by [3] [151 EE] 



p q {e) oc 



(4) 



with [x] + = x for x > and zero otherwise. For q — 1, Eq. (3) reduces to 

tf(e) = -e, (5) 

which leads to the Gaussian distribution given by 

Pi(e) oce~ (r/ ^ )e2 . (6) 

For g > 1 and g < 1, Eq. (4) yields long-tail and cut-off distributions, respectively. Thus 
Eqs. (2) and (3) generate the Gaussian and non-Gaussian noises, depending on the value 
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of parameter q. Expectation values of e and e 2 are given by 



(4 = 0, (7) 

which shows that (e 2 ) diverges at q = 5/3. 

In order to make our calculation tractable, we replace the e 2 term in the denominator 
of Eq. (3) by its expectation value: e 2 ~ (e 2 ), to get [H] 

K{e) ~ --, (9) 



by which Eq. (2) becomes 



re = -l—)e + <fn](t), (model A) (10) 



with 

r, = ^4- (ID 
(5-3?) ' 

A model given by Eq. (1) with Eq. (10) is hereafter referred to as the model A, which 
is discriminated from the model Ao given by Eqs. (l)-(3). The solid curve in Fig. 1 
expresses r q . We note that we get r q = 1 for q = 1, and r q < 1 (r q > 1) for q < 1 
(1 < q < 5/3). The dashed curve will be discussed in Sec. 3.1. 

Now we discuss the FPE of the distribution p(x,e,t) for Eqs. (1) and (10), which are 
regarded as the coupled Langevin model. We get 

^p(x,e,t) = --^-{[F(x) + e + I}p{x,e,t)} + '^-^p(x,e,t) 



1 d . . 1 /0\ 2 d 



e—ep{x,e,t) 



12) 



We define means, variances and covariances by 

(x m e n ) = J dx J dex m e n p(x,e } t). (m, n: integer) (13) 

By using the moment method for the coupled Langevin model [T3l [TJ], we get their 
equations of motion given by 

d(x) 



dt 



+ e + />, (14) 
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* r q r^ (15) 

f/V '"~ ; = 2(x[F(x) + e + I]) + ?p 2 , (16) 



dt 

d(e 2 ) 2 



21 A«vf*V. (it) 



(it r g r \ r 



d(xe) 
dt 



-- (e[F(x)+e + I}) -— (xe). 

r q r 


(18) 


and covariance denned by 




H = (x), 


(19) 


v = (e), 


(20) 


7 = (x 2 )-(xY, 


(21) 


C = (e 2 )-(e) 2 , 


(22) 


X = (xe)-(x)(e). 


(23) 



When we expand Eqs. (14)-(18) as x = \i + 5x and e = v + be around the mean values of 
fj, and u, and retaining up to their second order contributions such as ((5x) 2 ), equations 
of motion become 



dfi 
~dl 

dv v 



dt r q T 



/0 + /27 + " + (24) 

(25) 



dj 
~dt 



2(/ l7 + 0) + ^ 2 , (26) 



with 



dx e 

When we adopt the stationary values for u, ( and (j>: 



1/ ~ v s = 0, (30) 



equations of motion for /j, and 7 become 



2/17+ 



(33) 
(34) 



eft '" ' (1 - r,r/i) 

where r q is given by Eq. (11). It is noted that the stationary value of ( s = (2 — g)0 2 /r(5 — 
3q) in Eq. (31) is a little different from (e 2 ) = 2 /r(5 — 3q) in Eq. (8), which is due to 
an introduced approximation. 

We may express the effective DE for x as 



x 



F eff {x) + I eff (t) + a eff r]{t) + ipm, 



with 



feff{X) 
Ieff(t) 



F(x) 
lit), 



a 



eff 



1 - r q h 



(35) 

(36) 
(37) 
(38) 

(39) 
(40) 



from which Eqs. (33) and (34) are derived [131 E]- Equations (35) and (38) clearly 
express the effect of non-Gaussian colored noise. The effective magnitude of noise a e ff 
is increased with increasing q (Fig. 1). In contrast, with increasing r, the effective a e ff 
value is decreased for f\ < which is usually realized. 
The FPE of P(x,t) for Eq. (35) is expressed by 



= -| { (F«„ + /)P (l , t )} + I| 

i> 2 3 2 „, , 
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a eff —oc eff P{x,t) 



(41) 



which may be applicable to a e // depending on x [i.e. Eqs. (50) and (59)]. The stationary 
distribution is given by 



For Fix) 



lnP(x) = 2 
-Ax, we get 

P{x) qc exp 



dx 



eff 



+ 1 



^ 2 eff + ^\ 



2 \ 2 



A 



2 /(l + Ar 9 )+# 



x 



A 



(42) 



(43) 



For F(x) = ax — bx 3 , we get 
P(x) oc exp 



2 /(1 _ Tg(a _ 3M)+v , 2] 



' 2 bx 4 ^ 
ax h 2Ix 



(44) 



2.2 Comparison with other methods 



We will compare the result of the moment method with those of several analytical meth- 
ods: the universal colored noise approximation (UCNA) and functional-integral methods 
(FI-1 and FI-2). 
(a) UCNA 

The universal colored noise approximation (UCNA) was proposed by Jung and Hanggi 
[SI E] by interpolating between the two limits of r = and r = oo of colored noise, and 
it has been widely adopted for a study of effects of Gaussian and non-Gaussian colored 
noises. By employing the UCNA, we may derive the effective DE for the variable x. 
Taking the time derivative of Eq. (1) with ip = 0, and using Eq. (10) for e, we get 



x = F'x + e + 1, 

= ( F .-L\ i+ (Z±I 



(45) 
(46) 



When we neglect the x term after the UCNA, we get the effective DE for x given by 

x = F eff {x)+I eff {t) + a effV {t), (47) 

with 



F[x) 



a 



eff 



(1- 


- r q F>) 


(I 


+ rj) 


(1- 


~ r q F>) 




r q (j) 


(1- 


- r q F>) 



(48) 

(49) 
(50) 



where F = F(x), F' = F'(x), and r q and r q are given by Eqs. (40) and (11), respectively. 
It is noted that given by Eq. (50) generally depends on x, yielding the multiplicative 
noise in Eq. (47). 

For F(x) = —Ax, the stationary distribution is given by 

21 



P (x) oc exp 



A(l + Xr q ) / I 
& V A 



(51) 
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which agrees with the result of Eq. (43) with if) = 0. 
For F(x) = ax — bx 3 , we get 

a(l — ar q )x 2 6(1 — ^ar q )x A b 2 r q x 6 



P (x) oc [1 — T q (a — 3bx )] exp 



x exp 



2I e (l - aT q )x 2I c bT q x 2 

in in 



(52) 



whose functional form is rather different from that given by Eq. (44). 
(b) Functional-integral method (FI-1) 

Wu, Luo and Zhu [9] started from the formally exact expression for P(x, t) of Eqs. (1) 
and (10) with I(t) = given by 

J^OM) = -j£[F{x)P{x,t)] - ^{e(t)5(x(t)-x))-if;^{mH<t)-x)),(53) 

where (■) denotes the average over the probability P(x,t) to be determined. Employing 
the Novikov theorem [T7] and the functional-integral method, they obtained the effective 
FPE for P(x, t) which yields Eq. (47) but with 



F(x) 



a 



w 

eff 



r q <j) 



1 - r F' 
1 ' q 1 s 



(54) 
(55) 



where r q and r q are given by Eqs. (40) and (11), respectively, F' = dF/dx and F s et al. 
denote the steady-state values at x = x s . 
For F(x) = —Ax, we get 



P w (x) oc exp 
For F(x) = ax — bx 3 , we get 

P w (x) oc exp 



A(l + Ar g ) j2 



[1 - r q (a - 36/x 2 )] 



ax 2 — 



bx 4 



(56) 



(57) 



(c) Functional-integral method (FI-2) 

By applying the alternative functional-integral method to the FPE for p(x, e, t) given 
by Eqs. (1) and (10) with ip = I(t) = 0, Fuentes, Toral and Wio [1] derived the FPE of 
P(x,t), which leads to the effective DE given by Eq. (47), but with 

F 



fF 



eff 



a 



» 



eff 



;i - s q rF> 



;i - s q rF') 



(58) 
(59) 
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with 

(60) 

We note that afj-j generally depends on x, yielding the multiplicative noise in Eq. (47). 
For F(x) = —\x, we get 

P F (x) oc (1 + XTsJs-Wto-W exp 

with 



(62) 

For F(x) = ax — bx 3 , it is necessary to numerically evaluate the distribution P(x) with 
the use of Eqs. (42) and (58)-(60). 

A comparison among various methods is summarized in the Table 1. We note that 
the result of our moment method agrees with that of FI-1, but disagrees with those of 
UCNA and FI-2. The result of UCNA is not identical with that of FI-2, although they 
are consistent each other if the identity of s q = r q holds, which is realized for q — 1 with 

si = n = l. 

2.3 Numerical calculations 

We present some numerical examples to make a comparison with direct simulation (DS), 
which has been performed for Eqs. (l)-(3) by the fourth-order Runge-Kutta method 
with a time step of 0.01 for 1000 trials. Figures 2(a)-2(f) show the stationary probability 
calculated by various methods for F(x) = —Xx with changing parameters of q and r for 
fixed = 0.5 and tp — 0. A comparison between Fig. 2(c) and 2(d) shows that the width 
of the distribution for r = 1.0 is narrower than that for r = 0.5. This is explained by 
the reduced effective strength of ct e ff = 0/(1 + At) by an increased r. We note that for 
q = 1.0, results of all methods are in good agreement each other. Comparing Fig. 2(a) 
to Fig. 2(c) [and Fig. 2(b) to Fig. 2(d)], we note that the width of the distribution for 
q = 0.8 is a little narrower than that for q = 1.0. This is due to the fact that the r q 
value is reduced to 0.82 from unity. An agreement among various methods is good for 
q = 0.8. In contrast, Figs. 2(e) and 2(f) show that for q = 1.5, the width of P(x) becomes 
wider because of the increased ri. 5 = 2.0. The results of the moment method, UCNA and 
FI-1 are in fairly good agreement. On the contrary, the distribution calculated by FI-2 is 
sharper than that of DS. 



l + (<7 




Xr(q- l)s q 



(61) 



1 + (?-!) 



VA 2 



x 
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Figures 3(a)-3(f) show the stationary probability calculated by various methods for 
F(x) = x — x 3 with changing parameters of q and r for fixed = 0.5 and ip = 0.0. The 
general trend realized in Figs. 3(a)-3(f) is the same as in Figs.2(a)-2(f). The result of FI-2 
for q = 1.5 is not so bad compared to those of other approximation methods. However, 
the result of FI-2 for q = 0.8 and r = 1.0 is worse than other methods. 



3 Model B 



3.1 Moment method 

In order to generate non-Gaussian noise, we may employ an alternative model (referred 
to as the model B) given by 



x 
re 



F{x) + e{t)+I{t), 



(model B) 



(63) 
(64) 



where F(x) expresses an arbitrary function of x, I(t) an external input, r the characteristic 
time of colored noise, and a and (3 denote magnitudes of additive and multiplicative noises, 
respectively, given by zero-mean white noises, i] and £, with correlations: (r)(t)r)(t')) = 

(mm) = S(t - t') and (v(t)m) = o. 

The FPE for the distribution p(e, t) for Eq. (64) in the Stratonovich representation is 
given by 

9 , s 1 9 r , N1 1 fa\ 2 d ( d . . .A 1 f(3\ 2 d 2 , . 

= -g- e [^t)) + -(-) _( c _[q,( ei t)])+ 5 (5) ^P(e,t). (65) 



dt 



2 \tJ de \ de' 



2\r 3e 2 



The stationary distribution of e has been extensively discussed [I0]-[T1] in the context of 
the nonextensive statistics [T51 [16] - It is given by [T0]-[T4] 



p,(e) oc 



'OC 



-(r/a 2 +l/2) 



K(3 2 



9-1 



(66) 
(67) 



with 



1 + 



2a 2 



K 



2r + a 2 j 
3-q\ ( 2t 



2r + a 2 J ' 



(68) 
(69) 
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where [x] + = x for x > and zero otherwise. In the limit of a = 0.0 (q = 1), the 
distribution given by Eq. (67) reduces to the Gaussian distribution given by 



p(c) oc exp(--e 2 J. (70) 

In the opposite limit of (3 — 0.0, Eq. (67) leads to the power-law distribution given by 

p(e) oc e-\ (71) 

with 

5 = 1 + % = —. (72) 



cr q — 1 



The expectation values of e and e 2 are given by 



0. 



n(3 2 



ff 2 



(73) 
(74) 



r(5-3g) 2(r-a 2 )' 

The second moment is finite for a 2 < A (g < 5/3). It is expected that Eq. (64) leads to 
the non-Gaussian colored noise with the correlation given by 



«t)e(t')) 



P 2 



■ exp 



t-t' 



(75) 



2(r-« 2 ) 

By applying the moment method to Eqs. (63) and (64), we may obtain the effective 
one-variable DE for x given by 



* = Feff + I(t)+(3 e f f m, 



with 



F eff 
Peff 

U„ 



F(x), 



(3u q , 



1 



3-q 



(76) 

(77) 
(78) 
(79) 

(80) 



V(l-a 2 /r) V 5 - 3 ?' 
details of calculations being explained in the Appendix. The q dependence of u q is plotted 
by the dashed curve in Fig. 1, where u q < 1, u q — 1 and u q > 1 for q < 1, q — 1 and 
1 < q < 5/3, respectively. We note that u q has a similar g dependence as r q shown by the 
solid curve. 
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The FPE of P(x, t) for Eq. (76) is given by 







P{x,t) 







dt y ' ' dt 
The stationary distribution is given by 



[(F eff + I)P(x,t)} + 



d_ 

dx 



d 



(81) 



lnP(x) 
For F(x) = —Xx, we get 

P(x) oc exp 
For F(x) = ax — bx 3 , we get 

P(x) oc exp 



dx 



eff 



A 2 



In 



'A 2 



eff 



(82) 



eff 



'X(l + XrY 

, A a 2 , 



[1 - r(a - 36/i 2 )] 



ax 



b -^ + 2Ix 



(83) 



M) 



Equations (83) and (84) are similar to Eqs. (43) and (44) (with ijj = 0.0), respectively, for 
the model A, although u q and r in the former are different from r q and r q in the latter. 

It would be interesting to compare the result of the moment method for the model B 
with those of the UCNA and FI method, as we have made for the model A in Sec. 2.2. 
Unfortunately the UCNA method cannot be applied to the model B because Eq. (64) 
includes the multiplicative noise [TB|. It is very difficult to apply the FI method to the 
model B including both additive and multiplicative noises: such calculations have not 
been reported as far as we are aware of. Then we will make a comparison of the result of 
the moment method only with that of DS in the next subsection 3.2. 

3.2 Numerical calculations 

We present some numerical examples to make a comparison with DS, which has been 
performed for Eqs. (63) and (64) by the Heun method with a time step of 0.001 for 1000 
trials. Figures 4(a)-4(d) show the stationary probability P(x) calculated for F(x) = —Xx 
with changing parameters of a and r for a fixed (3 = 0.5. In Figs. 4(a) and 4(b) for 
a = 0.0 (q = 1.0), we observe that the width of P(x) is decreased with increasing r, 
as shown in Figs. 2(c) and 2(d). Figures 4(c) and 4(d) show that with increasing a to 
0.5, we get wider width in P(x) because we get q = 1.40 and 1.22 for r = 0.5 and 1.0, 
respectively [Eq. (68)]. 

Similarly, Figs. 5(a)-5(d) show P(x) for F(x) = x — x 3 . Results of the moment method 
are in good agreement with those of DS for a = 0.0 (q = 1.0) as shown in Figs. 5(a) and 
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5(b). The width of P(x) in Fig. 5(c) for a = 
that in Fig. 5(a) for a = 0.0 and r = 0.5 (q = 
5(d) for a = 0.5 and r = 1.0 (q = 1.22). 



0.5 and r = 0.5 (q = 1.40) is wider than 
1.0), but it is narrower than that in Fig. 



4 Discussion 



We will make a comparison among the various methods for generating non-Gaussian noise 
given by 



with 



x = F(x) + e(t) + I(t), 



(85) 



re = K(e) + (f)r)(t), (model A ) (86) 

re = -(^-J +(f)r}(t), (model A) (87) 

re = -e + e ar](t) + /3£(t), (model B) (88) 

where 77 and £ are white noises, r q is given by Eq. (11), and K(e) is given by Eq. (3) or 

dU(e) 



K(e) = 
U(e) = 



de 

At 



In 



l + (<7 



r 



(89) 
(90) 



2r(g - 1) 

Note that the model A is derived from the model A with the approximation: K(e) ~ 
— e/r q and U(e) ~ t 2 /2r q [Eq. (9)]. Noises in the models A and A are generated by a 
motion under the potentials given by Eq. (90) and U(e) = e 2 /2r q , respectively, subjected 
to additive noise. In contrast, noise in the model B is generated by a motion under the 



potential of U (e 



l /2 subjected to additive and multiplicative noises. 



We note from Eqs. (4) and (67) that the stationary distributions of e in the models 
A and B become the equivalent non-Gaussian distribution if the parameters in the two 
models satisfy the relation: 

P 2 



k(3 z 



for q > 1 



(91) 



(l + a 2 /2r)' 

This equivalence, however, does no hold between the models A and B, because the sta- 
tionary distribution of the model A is not the non-Gaussian but the Gaussian given by 



p q (e) oc exp 



r q 4> 2 



(model A) 



(92) 
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As for the dynamical properties, equations of motion for (e 2 ) in the moment method 
are given by Eqs. (17) and (A6): 

^ - -(^)«+g) 2 , (model A) ,93) 
M _ _| + (*)\ (modd B) 

Equations of motion for fj, and 7 are given by Eqs. (33), (34), (A16) and (A17): 

^ = / + / 27 + /(t) ) (models A and B) (95) 



dt 

— = 2/17 + 

(1 - r,r/i) 



37 = 2/17+ 7; 9 FT ' (model A) (96) 



where w„ = yJ(3-q)/(5-3q) [Eq. (A18)]. In the model A, we have adopted the approx- 
imation: if(e) ~ — e/r g [Eq. (9)], without which reasonable results are not obtainable in 
the moment approach (see the discussion below). Equations (93)-(97) show that equations 
of motion for the models A and B have the same structure. In the case of weak noise and 
small r, for which the second-moment approach is expected to be valid, the dynamical 
properties of the models A and B (as well as the model A ) are qualitatively the same, 
although there are some quantitative difference among them: i.e. the stationary value of 
7 of the model A is different from that of the model B. 

Our discussion presented in this paper is based on the second-order moment method. 
Effects of higher-order moment neglected in our method are examined in the following. 
The equation of motion for the fcth moment with even k of the model A is formally given 
by 



dt 



*) (e^K{e)) + k -^ Y ^ (^j (e^), (model A ) (98) 



though an evaluation of the first term of Eq. (98) is very difficult. In order to get 
a meaningful result within the moment method, we have assumed K(e) ~ —^l r q (the 
model A) to get 

' ( - I (c (A - 2) ). (model A) (99) 



d ^ k ) _ ( k \ , k(k-l) (o\ 2n 



dt \ r 'q T I 2 \T 
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from which we may recurrently calculate the stationary kth moment as 

= (J ^ i (^) <<*- 2) >, (100) 

_ (modelA) m 

The stationary distribution in the model A given by Eq. (4) leads to the second- and 
fourth-order moments: 

^ - < io2 > 

«*> - (s-frff-s^ - (modelAo) (103) 

In contrast, the stationary distribution in the model A given by Eq. (92) yields 

(f4> = 3<e '> 2= (5-"4)V' - (m ° delA) (105) 
The expression of Eq. (104) is different from that of Eq. (102) by a factor of (2 — q). 
The ratio of Eq. (105) to Eq. (103) becomes (2 — q) 2 (7 — 5q), which is less than unity 
for 1 < q < 5/3. These show that the distribution given by Eq. (92) in the model A 
underestimates the effective width of the distribution of e compared to that in the model 
A . In order to include the higher-order moment in an appropriate way, we have to go 
beyond the approximation with K(e) ~ — e/r q [Eq. (9)]. 

With the model B, we may obtain the equation of motion for the kth moment with 
even k, as given by 



d{e k 



2 , x 2" 1 ' - ' , . . , \ 2 



k k /a\ 



dt 

The stationary value of the fcth moment is given by 

(k - l)f3 2 jk-2) 



{t k ) + k{k 2 l) f^J (e<*- 2 >). (model B) (106) 



(e k ) = " (e^) (107) 

[ ' 2(r-ka 2 /2) [ h 1 } 



(k - 1)!! fi- 



ll 



2 fe / 2 nj£i(r - ia 2 

For example, second- and fourth-moments are given by 



(108) 



2[t — or J 

(e 4 ) = 77 T^C; 77, (model B) (110) 

x 1 4(r - 2a 2 )(r - a 2 )' ( J y ' 
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which agree with the result obtained from the stationary distribution given by Eq. (66) or 
(67). We get the positive definite (e k ) for a 2 < 2r/k. This suggests that for 2r/k < a 2 < r 
with k > 4, the kth moment diverges even if the second moment remains finite. This 
might throw some doubt on the validity of the second-moment approach. Equation (106) 
expresses that the motion of (e k ) depends on those of its lower moments (< k — 2), but 
it is independent of its higher moments (> k + 2). For example, even if (e 4 ) diverges, 
it has no effects on the motion of (e 2 ) for r/2 < a 2 < r. It is promising to take into 
account contributions from higher-order moments in the model B, although its validity 
range becomes narrower because a has to satisfy the condition: a 2 < 2r/k for the kth 
moment to remain finite. 

Our discussions presented in the preceding sections are confined to the stationary 
properties of the Langevin model subjected to non-Gaussian noise. It is possible to discuss 
its dynamics, by solving equations of motion for ji and 7. Numerical calculations for the 
model B are plotted in Figure 6(a) and 6(b), which show the time dependences of [i and 
7, respectively. We apply an external pulse input given by lit) — AQ(t — 100)0(200 — t) 
with A = 1.0, which is plotted at the bottom of Fig. 6(a), Q(x) denoting the Heaviside 
function. Figure 6(a) shows that n(t) of the moment method is in good agreement with 
the result of DS. Figure 6(b) shows that 7(t) is independent of an input pulse [Eq. (A17)]. 
With increasing a, a steady value of 7 is increased. The result of the moment method 
for a = 0.0 is in good agreement with that of DS although for a = 0.5, the former is 
underestimated compared to the latter. 

The overall behavior of the stationary distribution is fairly well reproduced by all 
the approximations mentioned in Sec. 2.2. Tails of P(x) are, however, not satisfactorily 
described, in particular, in calculations of the model A. This is partly due to the fact that 
the approximate Eq. (10) yields the Gaussian stationary distribution given by Eq. (92), 
though Eqs. (2) and (3) are originally introduced to generate non-Gaussian noise. This 
point is improved in the model B, in which the stationary distribution given by Eq. (64) 
is non-Gaussian as expressed by Eqs. (66) and (67). Indeed, tails of P(x) of Fig. 5 for 
the model B are slightly well reproduced than those of Fig. 3 for the model A. 

5 Conclusion 

To summarize, we have studied effects of non-Gaussian noise on the Langevin model, by 
using the second-order moment approach. The obtained result is summarized as follows. 
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(1) With increasing r, the width of the stationary distribution P(x) is decreased. 

(2) For q > 1 (g < 1), the width of P(x) is increased (decreased) compared to that for 
q=l. 

(3) The prefactor of F e ff for the model A in the moment method agrees with that in FI-1, 
but disagrees with that in the UCNA and FI-2 (Table 1). 

The items (1) and (2) are realized in both the models A and B. This may be explained 
by the q- and r-dependent a e ff [Eq. (38)] or j3 e jf [Eq. (78)]. As for the item (3), it 
is necessary to point out that although the UCNA [6j |7j exactly interpolates between 
the two limits of r = and r = oo, it is not exact for 0(t) [TH]. The functional 
integral method is a formally exact transformation if the functional integral is correctly 
performed. In the actual applications, however, it is inevitable to adopt some kinds of 
approximation, with which the final result depends on the adopted approximation. The 
difference between the results of FI-1 [5] and FI-2 [I] arise from the difference between the 
adopted approximations in performing the functional integral. These yield the difference 
in the results listed in the Table 1. 
As for the models A and B, we get 

(i) although the stationary distribution of p(e) in the model A is the Gaussian, the effect 
of the non-Gaussian distribution of the original model A is fairly well taken into account 
by a factor of r q , and 

(ii) the newly introduced model B, which yields the stationary non-Gaussian p(e) equiva- 
lent to that of the model A , is expected to be a promising model generating non-Gaussian 
noise. 

It is possible to apply the moment approach to a wide class of stochastic systems 
subjected to non-Gaussian noise, because its calculation is simple and transparent. It 
would be interesting to investigate effects of non-Gaussian noise on the synchronization 
in coupled nonlinear systems with the use of the model B, which is left as our future 
study. 
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Appendix 

We discuss an application of the moment method to the model B given by Eqs. (63) 
and (64), for which the FPE of the distribution p(x, e, t) in the Stratonovich representation 
is given by 



d d Id 

—p(x,e,t) = -—{[F(x)+e + I]p(x,e,t)} + - — [ep(x,e,t)\ 



1 fa\ 2 d 
+ 2(7 de 



9 , 
e—ep{x,e,t) 



1 (PY d 2 . 

+ 2 7 or 2P ^ e ^ 



(Al) 



We define means, variances and covariances by 



(x m e n ) = J dx J de x m e n p(x, e, t). (m,n: integer) 



(A2) 



After simple calculations using Eqs. (Al) and (A2), we get their equations of motion 
given by [131 EE 



d(x) 
dt 

<m 

dt 
d(x 2 ) 

dt 
d(e 2 } 

dt 
d(xe) 

dt 



(F(x) + e + I), 

2(x[F(x) + e + I\), 
2 /-^ 2 



T 



(e[F(x) + e + I\) -~(xe). 



(A3) 
(A4) 
(A5) 

(A6) 

(A7) 



We will consider the variables of fi, u, 7, ( and defined by Eqs. (19)- (23). Their 
equations of motion become 



dfi 

~dt 
dv 

~dt 

o?7 

~dt 

d( 

dt 
d<j> 
~dt 



/0 + /27 + " + 

v 1 (a\ 2 
-r + 2 It) 



2(/i7 + 0) 
2 . „ (a 



T \T J \T 



a 



T 



iw fa 



(A8) 
(A9) 
(A10) 

(All) 

(A12) 



where f t = {1/ V.)d l F{n) / dx £ . 



18 



When we assume the stationary values for u, ( and 0: 



C - C, = ^-^, (M4) 



v ~ v s = 0, (A13) 
2(r-a\ 

(j) - <Ps = 7w3 ¥r -T7Z TT, (A15) 

2(1 -a 2 /r)(l -r/i) 
equations of motion for /i and 7 become 

^ = /0 + /27 + /W, (A16) 

with 



«■ - V(T3P7o=fe (A18) 

With increasing r, m 9 is decreased because of a decreased g (Fig. 1). Equations (A17) 
and (A18) lead to the stationary value of 7 S given by 

l 2 (3~g) _ r , 

% " (-2A)(l-rA)(5-3 ? ) (-/OCl-r/O^ 1 J 

It is noted that equations of motion given by Eqs. (A16) and (A17) may be derived from 
the one- variable DE given by Eq. (76) [T3J [H]. 
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F eff 


J eff 


a eff 


method 


F 
F 

F/(l - r q rF') 
F/(l - s q rF>) 


I 

(I + ri)/(l-r q TF>) 


r q <Pl(\l l - r q r fl) 

r q <f>/Ql - r q rF' s ) 
r q (j)/{l - r q rF') 
Sq( j>/(1 - s q rF') 


moment ^ 

FI-1 2 ) 

UCNA 3) 
FI-2 4 ) 



Table 1 A comparison among various approaches to the model A [Eqs. (1) and (10)] 
yielding the effective differential equation given by x = F e ff + I e ff + a e // 77 (i), where 
r q = 2(2 - g)/(5 - 3q) and s q = 1 + (g - l)(r/20 2 )F 2 ; (1) the moment method: (2) 
functional-integral (FI-1) method of Ref. [9J: (3) UCNA calculation after Ref. jSHZj: (4) 
functional- integral (FI-2) method of Ref. [1] (see text). 
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Figure 1: The q dependence of r q [Eq. (11): solid curve] and u q [Eq. (80): dashed curve]. 



Figure 2: The stationary probability P(x) for F(x) = —\x of the model A [Eqs. (1)- 
(3)] calculated by direct simulation (DS: dashed curves), and P(x) of the model A [Eqs. 
(1) and (10)] calculated by the moment method (solid curves), UCNA (chain curves) 
and FI-2 (dotted cures): results of FI-1 agree with those of the moment method: (a) 
(q, t) = (0.8, 0.5), (b) (0.8, 1.0), (c) (1.0, 0.5), (d) (1.0, 1.0), (e) (1.5, 0.5), and (f) (1.5, 1.0) 
(0 = 0.5 and ip = 0.0). 



Figure 3: The stationary probability P(x) for F(x) = x — x 3 of the model A [Eqs. (1)- 
(3)] calculated by direct simulation (DS: dashed curves), and P(x) of the model A [Eqs. 
(1) and (10)] calculated by the moment method (solid curves), UCNA (chain curves) 
and FI-2 (dotted cures): results of FI-1 agree with those of the moment method: (a) 
(q, r) = (0.8, 0.5), (b) (0.8, 1.0), (c) (1.0, 0.5), (d) (1.0, 1.0), (e) (1.5, 0.5), and (f) (1.5, 1.0) 
(0 = 0.5 and ^ = 0.0). 



Figure 4: The stationary probability P(x) for F(x) = —Xx of the model B given by Eqs. 
(63) and (64) with (a) (a,r) = (0.0,0.5) (q = 1.40), (b) (0.0,1.0) (q = 1.22), (c) (0.5,0.5) 
(q = 1.40), and (d) (0.5, 1.0) (q = 1.22), calculated for f3 = 0.5 by DS (dashed curves) 
and the moment method (solid curves). 



Figure 5: The stationary probability P(x) for F(x) = x — x 3 of the model B given by 
Eqs. (63) and (64) with (a) (a,r) = (0.0,0.5) (q = 1.40), (b) (0.0,1.0) (q = 1.22), (c) 
(0.5,0.5) (q = 1.40), and (d) (0.5, 1.0) (q = 1.22), calculated for (3 = 0.5 by DS (dashed 
curves) and the moment method (solid curves). 



Figure 6: The time dependence of (a) fi(t) and (b) ^(t) of the model B for a = 0.0 
and a = 0.5 with (5 = 0.5 and r = 1.0, calculated by the moment method (solid curves) 
and DS (dashed curves), an input pulse being plotted at the bottom of (a). Results for 
a = 0.5 in (a) is indistinguishable from those for a = 0.0. 
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